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ABSTRACT 


During  hypervelocity  flight,  the  heat  transfer  rate  to  a  space  vehicle  depends 
on  the  contributions  of  aerodynamic  (convective,  conductive,  diffusive  and  chemically 
reactive)  transport  phenomena  and  radiative  transport  processes.  For  certain  com¬ 
binations  of  vehicle  flight  speed,  altitude,  vehicle  geometry  and  planetary  gas  compo¬ 
sition,  it  is  possible  for  the  radiative  transport  to  be  of  comparable  magnitude  to,  or 
even  to  predominate  over,  the  aerodynamic  processes,  and  hence,  the  twin  problems 
of  aerodynamic  and  radiative  heat  transfer  should  be  treated  self-consistently. 

The  macroscopic  equations  for  the  conservation  of  mass,  momentum  and 
energy  ( including  radiative  contributions )  are  first  presented.  The  hypervelocity 
flow  field  is  then  considered  to  consist  of  three  definable  regions;  namely,  the  shock 
wave,  the  influscid  flow  and  boundary  layer. 

The  equations  are  treated  for  the  optically  thin  and  optically  thick  limits  and 
numerical  results  are  obtained  for  the  effects  of  radiation  on  the  structure  of  the 
shock  wave,  and  the  structure  of  the  boundary  layer,  including  some  effects  on  heat 
transfer  rate  to  the  vehicle. 

In  order  to  assess  the  interplay  between  radiative  energy  and  chemical  energy, 
the  computations  of  shock  wave  structure  in  an  optically  thin  radiating  gas  are  carried 
out  for  two  limiting  cases.  In  the  first  model,  the  chemical  kinetic  processes  are 
assumed  to  be  initiated  after  the  shock  transition,  and  in  the  second,  the  chemical 
kinetic  processes  are  assumed  to  be  so  rapid  that  the  gas  is  everywhere  in  a  state  of 


dissociation  equilibrium. 


1.  INTRODUCTION 


The  theoretical  treatment  of  radiation  in  high  temperature  gases  over 
the  full  spectral  range  was  put  on  a  sound  basis  at  the  turn  of  the  present  century 
with  the  appearance  of  Planck's  law  of  radiation,  see  Ref.  (1).  A  treatment  of  the 
general  theory  of  radiative  transport  has  been  given  in  numerous  books  on  astrophy¬ 
sics,  Refs.  (2)  to  (4).  Methods  of  solution  of  the  radiative  transport  equation  includ¬ 
ing  the  effects  of  scattering,  (scattering  is  not  of  much  importance  in  the  problems 
of  interest  to  us  here),  are  given  for  various  conditions  in  the  excellent  books  by 
Chandrasekhar,  Ref.  ( 5 )  and  KOurganoff ,  Ref.  (6). 

Recently,  high  temperature  phenomena  have  become  of  increasing 
interest  to  engineers  and  scientists  concerned  with  the  problems  of  hypersonic  ballis¬ 
tic  re-entry  and  superorbital  entry,  into  the  earth's  atmosphere,  as  well  as  into  the 
atmospheres  of  the  neighboring  planets.  It  has,  therefore,  become  necessary  to  adapt 
the  theory  developed  principally  by  astrophysicists  to  the  conditions  of  interest  in 
hypersonic  aerospace  flight  problems,  where  one  of  the  major  differences  is  the  introduc¬ 
tion  of  a  solid  boundary;  namely,  the  vehicle.  We  begin  by  noting  that  a  number  of 
survey  articles  have  recently  appeared,  which  serve  to  describe  the  quantitative  state 
of  our  understanding  of  hypersonic  aerodynamic  and  radiative  heat  transfer,  see 
Refs.  (7)  -  (9),  and  hence,  we  will  restrict  our  discussion  here  to  conceptual  consi¬ 
derations. 

The  subject  of  radiative  heating  in  high  temperature  flowing  gases  can 
be  classified  in  terms  of  the  extent  of  the  coupling  between  the  radiative  transfer  and 
other  energy  transport  processes,  such  as  the  macroscopic  motion  of  the  material 


particles  and  the  chemical  reactions  in  the  gas. 

In  the  earliest  theoretical  treatments  of  radiation  at  hypersonic  speeds, 
(see  Figure  1),  it  was  assumed  that  the  shock  wave  (Region  I)  could  be  treated 
as  a  discontinuity,  and  that  the  viscous  boundary  layer  (Region  m)  was  very  thin. 

Since  the  equilibrium  gas  temperatures  for  suborbital  ballistic  entry  into  the  earth's 
atmosphere  are  usually  less  than  7000°K,  it  was  also  assumed  that  the  inviscid  shock 
layer  (Region  II)  was  optically  thin,  and  that  the  effect  of  radiative  emission  on  the 
gas  temperature  could  be  neglected.  Accordingly,  the  problem  became  one  of 
identifying  the  radiative  processes  in  high  temperature  air  and  calculating  the  frequency 
and  temperature  dependence  of  the  absorption  coefficients,  utilizing  experimental  data 
to  supplement  the  theory  (e.  g. ,  Refs.  (10)  -  (15)  ).  This  same  approach  has  also  been 
recently  applied  to  the  prediction  of  radiative  heat  transfer  at  hypersonic  speeds  in 
the  atmospheres  of  Mars  and  Venus,  Ref.  (16),  and  at  higher  gas  temperatures  during 
entry  into  Earth's  atmosphere,  Refs.  (17),  (18).  The  effect  of  departures  from  chemi¬ 
cal  equilibrium  on  radiative  heating  has  also  been  treated  within  the  framework  of  the 
assumptions  of  radiative  uncoupling,  and  an  optically  thin  gas  layer.  Refs.  (19)  -  (21). 

Attempts  are  also  being  made  to  consider  the  coupling  between  radiative 
transport  and  the  flow  processes;  and  departures  from  the  assumption  of  an  optically 
thin  gas  have  been  introduced  in  treating  the  hypersonic  shock  layer.  Because  of 
the  complicated  nature  of  the  problem,  attention  is  being  given  to  small  disturbance 
theory,  so  that  the  equations  may  be  linearized,  (e.g. ,  Refs.  (22)  -  (24) ). 

The  coupling  of  radiation  and  convection  has  been  treated  for  the  one- 
dimensional  case  by  the  Goulards,  Ref.  (25)  for  a  radiating  layer  adjacent  to  a  surface. 


In  their  paper,  they  extend  the  general  equations  developed  by  astrophysicists  to  Include 
wall  effects,  but  in  their  calculations  make  the  optically  thin  approximation. 

The  coupling  between  radiation  and  convection  in  an  inviscid  shock  layer 
has  been  treated  by  Goulard,  Ref.  (26)  and  also  by  Yoshikawa  and  Chapman,  Ref.  (27), 
using  the  general  equations  for  a  gray  gas. 

The  coupling  between  radiation  and  convection  in  a  viscous  layer  has 
been  treated  for  the  normal  shock  wave  when  the  gas  is  optically  thick  by  Sen  and  Guess, 
Ref.  (28)  and  by  Marshak,  Ref.  (29).  A  study  of  the  radiation-resisted  shock  wave 
has  been  carried  out  by  Clarke,  Ref.  (30). 

In  view  of  the  experimental  data  concerning  precursor  radiation  from 
shocks  in  air,  Ref.  (31),  attempts  are  being  made  to  treat  the  problem  analytically. 

A  study  of  the  upstream  photoionization  produced  by  a  strong  shock  wave  has  been 
carried  out  recently  by  Ferrari  and  Clarke,  Ref.  (32). 

The  interaction  between  radiation  and  convection  in  a  boundary  layer 
has  been  treated  by  Koh  and  De  Silva,  Ref.  (33)  for  the  viscous  flow  over  a  flat  plate, 
when  the  gas  is  assumed  to  be  optically  thin.  Viskanta  and  Grosh,  Ref.  (34)  have 
treated  the  wedge  boundary  layer  with  coupled  radiation  in  the  optically  thick  approxi¬ 
mation.  However,  they  evidently  did  not  realize  that  it  is  the  Rosseland  mean  absorp¬ 
tion  coefficient  which  is  appropriate  in  their  case,  rather  than  the  Planck  mean. 

The  most  recent  studies  of  radiation  coupling  in  the  hypersonic  boundary 
layer  include  the  work  of  Howe,  Ref.  (35),  who  considered  the  injection  of  an  absorbing 
gas  from  the  surface  of  the  vehicle,  and  the  study  of  the  merged  viscous  radiating 
layer  by  Howe  and  Viegas,  Ref.  (36),  utilizing  the  general  equations  in  the  gray  gas 
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approximation, 


It  is  also  appropriate  to  mention  that  theoretical  studies  of  the  spectral 
absorption  coefficient  of  air  and  air-like  species  are  being  carried  out  without  special 
reference  to  space  vehicle  applications,  (e.g. ,  Refs.  (37),  (38)  and  (39)  ). 

In  the  present  paper,  the  authors  will  consider  both  shock  wave  structure 
and  boundary  layer  structure,  for  the  optically  thin  and  optically  thick  approximations. 
As  a  point  of  departure,  one  of  the  authors  has  recently  re-derived  the  equations  for 
radiative  transport  and  the  coupling  between  radiation  and  matter,  for  quite  general 
conditions.  Ref.  (40),  utilizing  a  somewhat  different  point  of  view  than  is  usually 
employed  in  order  to  clarify  the  nature  of  the  approximations  which  go  into  the  usual 
form  of  the  equations.  In  this  reference,  the  particular  forms  of  the  radiative  contri¬ 
bution  to  the  macroscopic  equations,  which  apply  for  the  conditions  of  interest  to 
aerospace  flight,  are  obtained  from  the  general  equations. 
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2.  SYMBOLS 


a,  b,  c,  d  constants 

3  2  r  V1 

By(T)  =  2hv  /c  1  -  exp  (hy/ftT)  ,  Planck  Intensity 

Jeo  q  4 

B  (T)dt'  =  —  T* 
o  V  IT 


c 

c 


C 


i 


C 

P 


LC.C 
i  1 


EC  C 
i  1 


e 

E 


TV 


Ei*  Ej 


f  = 

n 


grgj 


velocity  of  light 

dimensionless  specific  heat  at  constant  volume 
mass  fraction  of  species  i 

specific  heat  at  constant  pressure 

frozen  specific  heat  of  the  gas  at  constant  pressure 

frozen  specific  heat  of  the  gas  at  constant  volume 

binary  diffusion  coefficient 

thermal  diffusion  coefficient 

internal  energy 
radiative  energy  density 

radiative  energy  density  spectrum 

energies  of  states  i  and  j 
similarity  stream  function 

dimensionless  tangential  velocity 

degeneracies  of  states  i  and  j 
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1 

% 


h 

hi 

Ah! 


I  (s) 


k 

ft 

K 

K 

V 

K 

P 

K 

P 

Kr 


V 


PC  D.. 
P»J 


Vw 


m 

M. 


M  =  JXM 
i  i  1 


n 

N 


a 

N.,  N 


j 


•  -  ^iipi 
v 

5*3 

»  ^ 


Planck's  constant;  static  enthalpy  of  mixture 

static  enthalpy  of  species  i,  including  chemical  enthalpy 

standard  heat  of  formation  of  species  i  evaluated  at  T  ^ 

specific  radiative  intensity  in  direction  s 

frozen  thermal  conductivity  of  mixture 

Boltzmann  constant 

mean  absorption  coefficient 

absorption  coefficient 

Planck  mean  absorption  coefficient 

equilibrium  constant  of  reaction  i 

Rosseland  mean  absorption  coefficient 

thermal  diffusion  ratio 

frozen  Lewis  number 

dimensionless  product  of  density  and  viscosity 

integration  constant  in  continuity  equaticm 
molecular  weight  of  species  i 
mean  molecular  weight  of  mixture 
unit  vector  in  the  negative  y  direction 
density  of  absorbers  of  type  "a" 
occupation  numbers  for  energy  states  i  and  j 
static  pressure 
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I 


integration  constant  in  momentum  equation 


frozen  Prandtl  number 


radiation  path  length 

integration  constant  in  energy  equation 


heat  flux  vector 


radiative  energy  flux  vector 


vector  radiative  energy  flux  spectrum 


dimensionless  density  (shock  wave) 
universal  gas  constant 


gas  constant  of  species  i 
nose  radius  of  body 

unit  vector  in  the  direction  of  photon  propagation 


temperature 


x  component  of  velocity 
y  component  of  velocity 
macroscopic  stream  velocity 
free  stream  velocity 


chemical  source  term,  mass  rate  of  production  of 
species  i  by  chemical  reaction  per  unit  volume,  per 
unit  time 

mole  fraction  of  species  i 
dimensionless  distance 


constant  of  integration  in  species  equation 


I 


body-oriented  coordinate  system 


x*  r0 
a,  fi 

6 

e 

w 

e 

w  v 


V,  t 

e 

M 

P 

cj 

a 

a  <v) 
a 

T 

V) 

n 

Subscripts 

A 

e 

i 

M 

r 


T 


constants 

thickness  of  gas  layer 

total  (hemispherical)  emissivity  of  the  wall 

spectral  emissivity  of  the  wall 

similarity  variables 

cos  0  =  n  •  _s;  dimensionless  temperature  (boundary 
layer) 

viscosity  of  mixture 
density 

Stefan -Boltzmann  constant 

dimensionless  pressure  (shock  wave) 

cross  section  for  absorption  of  photons  of  frequency 
v  by  absorbers  of  type  "a" 

dimensionless  temperature  (shock  wave) 
dimensionless  velocity  (shock  wave) 
solid  angle 


atoms 

outer  edge  of  boundary  layer 
th 

i  species 
molecules 
radiation 
stagnation  point 
thermal 


wall 


upstream 

downstream 

denotes  differentiation  with  respect  to  rj 
frequency 

reference  state 

thermal 

dimensionless 
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3.  GOVERNING  EQUATIONS  FOR  RADIATING  FLOWS 

When  the  approximation  of  local  thermodynamic  equilibrium  is  made 
and  scattering  is  neglected,  the  radiative  transport  equation  becomes 


SI  (s)  |-  -i 

+  s  •  V  I  (s)  =  K  I  B  (T)  -  I  (s)  I  . 
C  S  t  -*  V  -»  V  I  v  v  -*  I 


(3. 1) 


In  this  equation,  K 


V 


is  the  absorption  coefficient  defined  by  the  equation 


K 

V 


EN  or  (y) 
a  a  a' 


[ 


1  -  exp  (-  hy/  *T) 


]• 


(3.2) 


Here,  N  is  the  number  of  absorbers  of  type  "a"  per  unit  volume  and  has  the  value 
Si 

characteristic  of  thermal  equilibrium  at  the  local  matter  temperature  T,  while  o  (y) 
is  the  cross  section  for  absorption  of  photons  of  frequency  y  by  these  absorbers.  The 
factor  -  exp  ( -  h  y  /  *T  )^J  in  Eq.  (3.2)  has  the  effect  of  including  induced  emis¬ 
sion  with  true  absorption.  I^(s)  is  the  specific  intensity  of  radiation  of  frequency  y . 
The  quantity  I^(s)dydO  is  the  amount  of  radiative  energy  with  frequencies  between 
y  and  v  +  d  v  traveling  within  the  element  of  solid  angle  dfl  about  the  direction  of 
unit  vector  s  and  crossing  unit  area  perpendicular  to  s  per  unit  time.  The  thermal 
equilibrium  value  for  I ^ (s )  is  the  Planck  intensity  B^(T)  given  by  the  equation 


B  (T)  = 
V 


2hv' 


|exp  (hy/ftT)  -lj 


-1 


(3.3) 


The  time  derivative  term  in  Eq,  (3. 1)  is  negligible  unless  the  tempera- 
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f 


tures  are  very  high  ( It  T  ~  kev. ) ,  and  there  is  fairly  rapid  variation  of  tempera¬ 
ture,  density,  etc.,  with  time.  Thus,  in  problems  of  interest  here,  we  can  neglect 

this  term  and  write  the  radiative  transport  equation  as  follows: 


V  I  (S)  =K  B  (T)  -  I  (s)  . 
V  -  V  L  V  v  '  — 'J 


The  macroscopic  radiative  transport  equation  is  obtained  by  multiplying  Eqs.  (3.1) 
or  (3.4)  by  di/dO  and  integrating  over  v  and  solid  angle  f) .  The  result  in  the  latter 


v.  9r  =4nfdVKv  [b^T)  -cEri//4ffJ 


Here,  E^dp  is  the  energy  density  of  radiation  with  frequencies  between  v  and  p  +  dp, 
while  Qr  is  the  total  radiative  energy  flux  vector.  These  are  given  by  the  equations 


E  =  r  E  dp;  E  =  c  f  I  (s)dfl 
r  J  TV  rv  «  v 


(3.6) 


=  I  Srv*v,  9rv  = I  ■  y»>«>  • 

o  n 


where  the  integration  over  O  is  over  all  solid  angle,  i.e. ,  over  4 it.  For  the  sake  of 
completeness,  we  have  also  included  the  quantities  (the  flux  spectrum)  and  Ef 
(the  total  radiative  energy  density). 

The  approximation  of  local  thermodynamic  equilibrium  used  in  writing 


I 


Eqs.  (3.1)  and  (3.4)  means  that  the  microscopic  matter  distribution  functions  can  be 
approximated  by  those  characteristic  of  thermal  equilibrium  at  the  local  temperature 
T  in  writing  the  radiative  transport  equation.  That  is,  for  local  thermodynamic  equili¬ 
brium,  the  occupations  of  each  pair  of  energy  states  are  related  by  a  Boltzmann  factor 
as  in  the  case  of  true  thermodynamic  equilibrium 


'5r-T-exp{-<vEJ>/*T} 


(3.8) 


except  that  now  T  is  a  function  of  position  rather  than  being  a  constant.  In  this  equa¬ 
tion,  N.,  E.  and  g.  are  the  occupation  numbers,  energies  and  degeneracies,  respectively, 
of  the  state  i.  Symbols  with  subscript  j  refer  in  a  corresponding  manner  to  level  j. 

The  criteria  for  the  validity  of  the  approximation  of  local  thermodynamic 
equilibrium  in  the  treatment  of  the  radiation  is  that  the  interaction  of  the  matter  with 
the  radiation,  i.e. ,  absorption  and  emission  of  photons,  gives  a  small  contribution  to 
the  right  hand  side  (the  collision  part)  of  the  Boltzmann  equations  for  the  material 
particles.  But,  when  this  is  the  case,  it  is  also  valid  to  solve  the  material  particle 
Boltzmann  equations  with  the  interaction  with  the  radiation  neglected  and  obtain  the 
usual  form  for  the  macroscopic  energy  transport  equation  applicable  when  radiation 
is  unimportant.  However,  this  macroscopic  energy  transport  equation  is  not  correct 
until  either  side  of  Eq.  (3.  5)  is  included  in  the  usual  expression  for  the  divergence  of 
the  energy  flux  vector.  This  subject  is  discussed  by  one  of  the  authors  in  much  greater 
detail  in  Ref.  (40). 

The  approximation  of  local  thermodynamic  equilibrium  is  usually  valid; 
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if 


however,  in  the  case  of  bound -bound  absorption,  i.  e. ,  molecular  band  and  atomic 
line  absorption,  the  absorption  coefficient  is  proportional  to  the  matter  density  p , 
while  the  rate  of  collisions  between  material  particles  is  proportional  to  p  ,  Thus, 
at  very  low  densities,  interaction  with  the  radiation  dominates  over,  or  is  comparable 
with,  interactions  between  material  particles  in  determining  bound  state  occupation 
numbers.  In  this  case,  the  approximation  of  local  thermodynamic  equilibrium  breaks 
down  as  far  as  bound  state  occupation  numbers  are  concerned.  When  this  occurs,  the 
gas  is  likely  to  be  optically  thin  because  the  density  is  low,  which  leads  to  some  simpli¬ 
fications. 

In  the  present  paper,  we  always  assume  local  thermodynamic  equili¬ 
brium  to  exist  except  that  in  some  cases  the  chemical  composition  is  not  assumed  to 
be  that  characteristic  of  thermodynamic  equilibrium.  That  is,  we  assume  that  Eq.  (3.  8) 
applies  for  the  relative  occupation  of  the  energy  levels  of  a  particular  species,  but 
the  concentration  of  species  is  not  the  thermal  equilibrium  value.  The  only  effect  on 

the  radiation  equations  is  that  N  in  Eq.  (3.2)  should  be  changed  accordingly. 

3. 

In  order  to  obtain  the  explicit  values  for  E  ,  E  ,  Q  and  Q 

rv  r  ru  r 

appearing  in  Eqs.  (3.  5)  through  (3. 7)  one  must  solve  the  radiative  transport  equation, 
Eq.  (3.4).  The  general  solution  to  this  equation  is 

y<i>=J  vqi)BJT<qi>] 

Q  L  -J  a* 

1 

q  <  3*  9> 

+  I„(q0)  exp-J*  Kv(qi)dqx  , 
q« 
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where  q  is  the  position  along  the  direct! cm  s 


8  •  V  = 


(3.10) 


%  *  ql  *  q 


(3.11) 


qQ  is  the  location  along  s  of  a  boundary.  The  shocked  region  in  front  of  a  blunt  object 
traveling  at  hypersonic  speeds  in  a  planetary  atmosphere  can  be  approximated  by  a 
plane  parallel  layer.  If  y  is  taken  as  the  preferred  direction  and  we  choose  unit  vector 
n  to  be  in  the  negative  y  direction,  we  can  make  the  replacement 


dq  =  -  dy  /  cos  8  ,  cos  8  =  s  •  n  . 


(3. 12) 


Unless  the  gas  layer  is  either  optically  thick  or  optically  thin,  Eqs.  (3.5),  (3.6),  (3.7) 
and  (3.9)  are  difficult  to  treat  numerically  even  when  the  gas  under  consideration  can 
be  assumed  to  be  a  plane  parallel  layer  and  the  gray  gas  approximation, 

Ky  =  K  is  made.  The  latter  approximation  is  poor  for  high  temperature  air, 
as  is  discussed  later. 

In  the  present  paper  we  make  the  assumption  that  the  shocked  gas 
layer  is  either  optically  thick,  or  optically  thin.  The  conditions  that  must  be  met  for 
a  gas  layer  of  thickness  5,  if  it  is  to  be  optically  thick  or  optically  thin  are: 


J  K  (y)dy>>l;  optically  thick  layer 
o 


(3. 13) 


I 


(3. 14) 


J*  K  (y)dy<<  1;  optically  thin  layer, 
o  v 


First,  we  consider  the  optically  thick  case. 

3.1  Optically  Thick  Layer  of  Gas 

When  condition  (3. 13)  is  satisfied  for  all  important  frequencies  and  the 
change  in  temperature  within  a  space  interval  1  is  small,  one  can  make  a  Taylor 
expansion  of  1^  (s)  about  B^(T).  The  result,  which  is  seen  by  inspection  to  be  a  solu¬ 
tion  of  Eq.  (3.4)  is 

I  (s)  =  B  (T)  -  ■—  s  •  V  B  ( T )  +  ——  s  •  7  S  •  ^  B  (T) 

U  -  V  K  -»  V  K  -*  IK  -*  v  J 
V  V  V 

(3. 15) 


;  optically  thick  gas. 


For  the  one-dimensional  problem,  as  seen  by  Eqs.  (3.10)  and  (3.12), 


s  •  7  =  "  cos  6  d  /  dy  . 


(3.16) 


With  the  use  of  Eqs.  (3.6),  (3.7)  and  (3.15)  we  obtain 


E  .-AIL  b(T)  =  -4— -  ,  E  =—  B  (T), 
r  c  '  '  c  rv  c  j/ 


(3.17) 


where 
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(3.19) 


Thus,  in  an  optically  thick  gas  the  radiative  contribution  to  the  total  energy  flux  is  like 
a  thermal  conduction  term  with  a  coefficient  of  conductivity 


k  -  16aT3 

eff  3K^ 


(3.22) 


3.2  Optically  Thin  Layer  of  Gas 

When  condition  (3.14)  is  satisfied  the  exponential  factors  in  Eq.  (3.9)  can 
be  replaced  with  unity  so  that 

I  (q)  -  W  BrT(V1  +1  (^  .  (3.23) 

%  L 

Condition  (3. 14)  implies  that  the  first  term  on  the  right  hand  side  of  this  equation  is 
negligible  relative  to  B^(T)  .  Thus,  if  the  gas  is  optically  thin  for  all  important  fre- 
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quencies,  Eq.  (3.5)  becomes 

V.  ^-4,J-*d„K  [b  (T)  -  -Jri  Wdn]' 

o  “  O  J 


(3.24) 


qQ  is,  of  course,  a  function  of  angle,  or  s .  For  the  case  of  a  planar  wall  at  tempera¬ 
ture  T  on  one  side  of  the  gas  and  non-emitting  material  on  the  other  side,  the  approxi- 
w 

mate  situation  for  the  shocked  region  of  interest  here. 


f  I  (q  )dO 
4 n  ^  f'o' 


C  B  (T  ) 
W v  V  W 


(l  -  e, 


w>j' Vy)B,[T<y,"U  (3-25) 


The  first  term  is  the  contribution  due  to  emission  by  the  wall  material  and  the  second 

term  the  contribution  due  to  reflection  by  the  wall  of  radiation  emitted  by  the  gas. 

From  condition  (3. 14)  we  see  that  the  latter  term  must  always  be  negligible  relative 

to  K  B  (T)  and  thus  cannot  contribute  significantly  to  Eq.  (3.24).  If  we  assume  that 
V  v 

the  wall  emits  as  a  gray  body,  i.  e. ,  e  =  €  ,  where  e  and  e  are  the  spectral 

wv  w  wv  w 

and  total  hemispherical  emissivities  of  the  wall  we  get 


V-  Sc  -  4.  f  dvKv  [Bp<T)  -  -f  -  Bi)(Tw)] 

=  40[Kp(T)T4--f^Kp(T,Tw)Tw4]  . 


in  which  is  the  Planck  mean  absorption  coefficient  given  by  the  equation 


Kp(T)  =  J*  Kp  Bp(T)  dp/ B(T)  . 
1  o 


(3.26) 


(3.27) 


18 


The  reason  we  wrote  K  (T,  T  ) in  the  second  term  in  the  final  form  of  the  right  hand 
p'  w' 

side  of  Eq.  ( 3. 26)  is  because  is  a  function  of  the  gas  temperature  T,  while  we  inte¬ 
grated  over  B^(Tw),  i.e. , 


K  (T,  T  ) = 
p'  w 


J*  VT>  v\)dv 

o _ 

B<TW) 


(3.28) 


Actually,  in  the  studies  carried  out  in  this  paper,  the  wall  temperature  was  chosen  to 
be  small  enough  that  the  second  term  on  the  right  hand  side  of  Eq.  ( 3. 26)  was  negligible. 
When  the  assumption  is  made  that  e  =  the  net  flux  at  the  wall  in  the  optically  thin 
case  is  g 

Q  =  2o  J  K  T4dy-2(l  -f  )fff  K  T4dy-(  ctT4 
rw  0  p  w  Jp  ww 


=  e  {20  J6  K  T4dy  -  crT  41 
wl  «  p  w  ( 


(3.29) 


In  the  first  form  of  the  right  hand  side,  the  first  term  gives  the  half  of  the  emission  of 
radiation  by  the  gas  which  moves  toward  the  wall,  the  second  term  accounts  for  reflec¬ 
tion  of  part  of  this  radiation  by  the  wall,  and  the  third  term  is  the  wall  emission. 
Actually  as  is  discussed  in  Ref.  (40)  Eqs.  (3.25)  and  (3.29)  are  not  quite  correct 
unless  the  wall  reflectivity  and  emissivity  with  respect  to  intensity  I^(s)  are  inde¬ 
pendent  of  s  •  n  . 

Comparing  Eqs.  (3.21)  and  (3.27)  we  see  that  the  mean  absorption  coef¬ 
ficients  and  K  entering  the  expressions  for  optically  thick  and  optically  thin  gases 
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iKJ 


are  different  unless  depends  very  slightly  on  v.  For  constant  ,  they  are  of 
course  equal,  otherwise  >  K^.  In  the  case  of  an  air-like  gas,  the  frequency 
dependence  of  is  sufficiently  great  that  is  often  one  or  two  orders  of  magnitude 
greater  than  ,  as  seen,  for  example,  from  the  results  of  Stewart  and  Pyatt,  Ref.  (37) 
for  nitrogen.  Thus,  in  the  intermediate  case  when  the  gas  is  neither  optically  thick  nor 
optically  thin,  it  is  probably  nearly  always  a  very  poor  approximation  to  use  one  single 
average  value  for  the  absorption  coefficient  for  the  entire  spectrum,  (i.  e. ,  the  gray 
gas  approximation).  For  example,  in  using  K  =  in  the  general  equations  for  the 
gray  gas  approximation,  if  =  5,  one  would  approach  the  optically  thick  solution 
for  which  is  inappropriate.  On  the  other  hand,  for  the  same  conditions  of  tempera¬ 
ture,  density,  etc. ,  if  one  used  K  =  K^,  which  is  appropriate  in  the  optically  thick  limit, 

the  value  of  K  would  be  small,  perhaps  as  0. 1,  then,  the  general  equations  for  the 
R 

gray  gas  approximation  would  approach  the  optically  thin  solution  for  which  K  is 

K 

inappropriate.  In  fact,  the  variation  of  with  frequency  is  usually  large  enough, 
especially  at  the  higher  temperatures  of  interest  here,  that  unless  J  K^dy  is  very 
much  less  than  unity  or  J  IC^dy  very  much  greater  than  unity,  where  the  integration 
is  over  the  gas  thickness,  the  gas  is  likely  to  be  optically  thick  for  some  important 
spectral  regions  and  optically  thin  for  others .  Then  difficult  as  this  may  be,  in  order 
to  obtain  results  which  are  at  all  accurate,  it  is  probably  necessary  to  devise  differ¬ 
ent  appropriate  approximation  schemes  for  the  differing  spectral  regions. 

3. 3  The  Macroscopic  Equations 

Following  Hirschfelder,  Curtiss  and  Bird,  Ref.  (41),  the  governing 
macroscopic  equations  for  a  multicomponent  chemically  reacting  gas  may  be  written 
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In  the  following  form.  The  conservation  of  mass  is  given  by 


I 


It  is  appropriate  to  point  out  that  in  writing  the  equations  which  are 
applicable  in  the  optically  thick  limit,  i.  e. ,  when  the  condition  given  in  Eq.  (3. 13) 
is  satisfied  for  all  important  frequencies,  we  have  implicitly  assumed  that  the 
additional  condition 


VT  <<  1 


(3. 37) 


required  for  the  validity  of  the  expansion  (3. 15)  is  satisfied.  In  the  special  case 


that  the  region  under  consideration  consists  of  the  space  interval  Ay  =  y„  -  y  , 

Ct  i. 

between  the  position  y  at  which  the  temperature  Tj  is  very  low,  and  the  position 
y2  at  which  the  temperature  reaches  a  high  value  Tg,  as  is  usually  the  situation  in 
a  boundary  layer  or  shock  wave,  either  condition  (3. 13)  or  (3.  37)  implies  that  the 
other  condition  holds  approximately,  as  well.  This  is  seen  by  considering  Eq.  (3. 37) 


7T  L 


4(VT1> 

K  TAy 

V 


K  Ay 
v 


(3.  38) 


which  implies  condition  (3. 13).  Whether  or  not  condition  (3. 37)  Indeed  applies 
should  be  tested  after  the  numerical  solution  for  a  given  situation  is  obtained. 

One  might  expect  that  it  would  not  apply  for  the  shock  wave  (Region  I)  and  boundary 
layer  (Region  HI)  which  are  usually  very  thin  physically.  When  these  regions  are 
not  very  thin,  the  gas  density  is  usually  quite  low,  and  hence  it  is  expected  that  they 
would  still  be  optically  thin. 

In  the  optically  thin  limit,  Eq.  (3. 35)  applies  if  the  single  condition 
given  by  Eq.  (3. 14)  is  satisfied  for  all  important  frequencies.  However,  in  the 
latter  equation,  6  now  represents  the  total  thickness  of  the  shock  layer,  not  just 
the  extent  of  the  boundary  layer  or  the  shock  wave. 
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4.  RADIATION  IN  A  SHOCK  WAVE 

We  begin  by  noting  that  one  of  the  authors.  Refs.  (42),  (43)  has  studied 
the  structure  of  a  normal  shock  wave,  under  conditions  where  the  important  physico¬ 
chemical  coupling  was  the  excitation  of  the  internal  degrees  of  freedom  including  rota¬ 
tional  and  vibrational  modes.  The  theoretical  development  presented  here  follows  this 
earlier  work  in  some  features,  except  that  in  order  to  include  the  radiative  phenomena 
within  the  framework  of  a  relatively  simple  model,  internal  relaxation  coupling  was 
not  treated  at  this  time. 

4. 1  Optically  Thick  Shock  Wave 

Upon  utilizing  the  equations  presented  in  Section  3.  3,  the  governing 
equations  for  the  one-dimensional  steady  flow  of  a  viscous,  conducting,  diffusing, 
radiating  gas  are  obtained  as  follows.  The  conservation  of  mass  is  given  by: 


f  (PV)=0 

which  upon  integration  becomes: 


pv  =  m  (a  constant) 


The  conservation  of  momentum,  neglecting  radiation  pressure  is  given  by: 


du  _  dp  +  d  I  4  dv  1 
dy  dy  dy  3  ^  dy  J 


(4.3) 


which  upon  integration  becomes 


4  c3v 

mv  +  p  -  — —  u  -r~  =  P  (a  constant) 
3  dy 


(4.4) 


The  conservation  of  species  i  is  given  by: 


I 


! 


— -  (pvC,  +  j.)  =  w. 
dy  '  i  1  i 


which  when  integrated  becomes: 


mCj  +  ]^  -  4>.  =  Z  (a  constant) 


(4.5) 


where 


d$.  =  w.dy  , 
x  1 


(4.6) 


(4.7) 


and  the  one -dimensional  diffusion  flux  vector  for  a  binary  mixture  may  be  written, 
Ref.  (44): 


ii-P 


-  «  dCi  r(M-M.) 

dy  -2  L  M.  i  dy  lnP 


,  (  dC.  MM.  f"  i 

\— i-  +  — i-J-  - 

ij  l  dy  m2  L 

irl"T]| 


(4.8) 


+  K„ 


where  the  thermal  diffusion  ratio  is  given  by,  Ref.  (41) 

T 


D. 


K  = 
T 


n2M  M. 
i  1 


Jx. 

i] 


(4.9) 


For  an  optically  thick  gas,  the  energy  equation  may  be  derived  utilizing  Eqs. 
(3.  32),  (3.33)  and  (3.34): 


pv 


dh  d  f 

dy  dy 


,  ..  .  16CTT  . 
k(1  3kK„  * 


dT  ~1 

dy  J 


T“  (?jA)  + 

dy  i  i  i 


4  /  dv  \ 

T  **  VdT/ 


+  v 


j£- 

dy 


(4.10) 


Upon  introducing  the  momentum  equation,  the  definition  of  the  specific  heat  at  constant 
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volume  (assumed  constant)  and  then  Integrating,  we  obtain  for  the  conservation  of 
energy  in  a  dissociating  diatomic  gas: 


2mC  T  +  2pv  +  mv2  -  2k  (1  +  —  )  ~- 

v  r  '  3k  K  '  dy 


t--t  *»=. 


K  • 


<cv  -cv  )T 
A  VM  A 


!)a  [«' 


PK„ 


t2i‘ 


•] 

r] 


(4. 11) 


The  equation  of  state  is  given  by: 
p  =  pRT 

which  may  also  be  written: 


p  =  prm<1  +  ca)t 


(4. 12) 


(4.13) 


4. 2  Optically  Thin  Shock  Wave 

For  an  optically  thin  shock  wave  one  utilizes  Eqs.  (4,2),  (4.4),  (4.6) 
and  (4. 13)  for  the  conservation  of  mass,  momentum  and  chemical  species,  and  the  equa¬ 
tion  of  state  respectively,  but  now  one  utilizes  instead  the  following  form  of  the  conser¬ 
vation  of  energy  which  is  derived  from  Eqs.  (3.  32)  and  (3.  35): 


pv 


^  -  d  (k-^)-4K^oT4--^r  (Ejjh.  ) 


dy  dy 


dy 


Mir) 


+  v 


d£_ 

dy 


(4. 14) 


which  can  be  integrated  once  to  yield: 
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2  /IT  j  A 

2mC  T  +  2pv  +  m  -  2k  +  8  a  \  K  T  dy 

v  dy  “  p 


-f  *"-|-t2i”cAK0->cvM  .. 


!)a  [« 


PK_ 

t2i*  'VVT+i\0t^^A 


[A\°  -  (CV..  -  cv4> T] 

r]  ■ 


(4. 15) 


Q 


The  final  set  of  equations  for  both  the  optically  thick  and  optically  thin 
shock  wave  have  been  transformed  into  dimensionless  variables  in  Appendix  B. 

With  regard  to  the  chemistry,  we  note  the  following.  If  one  represents 
the  dissociation  equilibrium  of  an  air-like  molecule  by  the  reaction: 

M«*r2A  (4.16) 

then,  the  equilibrium  constant  can  be  written  in  terms  of  the  partial  pressures  of  the 
atoms  and  molecules,  or  in  terms  of  the  mole  fractions. 


K 


<PA>2  <xa)2p 


(4. 17) 


However,  since  the  mass  fraction  CA  appears  explicitly  in  the  governing  equations,  it 
is  convenient  to  introduce  the  relationship  between  the  mole  fraction  and  the  mass  fraction, 


and  then  it  can  readily  be  shown  that: 

K  1/2 

r 

c. 


„  K 

-  r _ 24 _ 1 

L  4P  +  K  J 


(4. 18) 


In  those  computations  where  the  gas  was  assumed  to  be  in  dissociation  equilibrium, 


the  equilibrium  constant  utilized  was  taken  as 


L*-'<!fc1»l>tt99ag,vat)iUg  "gw*? 


log 


10 


7.00 


9  x  104 
T  (OR) 


(4.19) 


which  corresponds  to  the  dissociation  of  nitrogen.  The  other  thermodynamic  and  trans¬ 
port  properties  of  nitrogen  were  based  on  data  appearing  in  Ref.  49  and  the  dimensionless 
groups  appearing  in  the  equations  were  calculated  in  a  manner  analogous  to  Refs.  42  and  43. 


5.  RADIATION  IN  A  BOUNDARY  LAYER 


In  this  section,  we  consider  the  effect  of  the  coupling  between  the 
radiation  and  the  viscous  flow  in  the  hypersonic  laminar  boundary  layer. 

5. 1  Optically  Thick  Boundary  Layer 

Upon  utilizing  the  equations  presented  in  Section  3.  3,  and  introducing 
the  Prandtl  boundary  layer  approximation,  one  obtains  the  following  governing  equa¬ 
tions  for  the  hypersonic  laminary  boundary  layer  of  a  viscous,  conducting,  diffusing, 
chemically  reacting,  radiating  gas  for  steady  flow  over  an  axi-symmetric  body.  The 
conservation  of  mass  is  given  by: 

~ir  (p>bu)+  <prov> =  0  (5-1) 


The  conservation  of  momentum  is: 
du  .  _  3u 


pu  IT  +  pv  ay 


ie. 

dx 


a  au 
ay  p  ay 


(5.2) 


where  it  is  noted  that  =  0  .  The  conservation  of  species  i,  assuming  the  appli- 

ay 

cability  of  Fick's  law  of  diffusion,  is  given  by: 


ac. 


pu 


i 


ac. 


+  pv 


a 

Sy 


fp  &  -ir-1  +  w  (5.3) 

L  i  j  ay  T  ay  J  i 


ax  ay 

Finally,  utilizing  Eqs.  (3.32),  (3.33)  and  (3.  34)  the  conservation  of  energy  may  be  written: 
P< 

j,  x  -j 

„3 


/  aTx  ST  \  _  iEL  ,  „  /_au_\ 
Cp  V"17*v17J-,,'Sr  M  V  Sy  / 


r  ..  .  160T"  .  JT 

+  dy  L  (  3kKn  )  ay 


■] 


(5.4) 
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j  fl  M2 
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which  applies  only  when  the  boundary  layer  is  optically  thick. 

Upon  introducing  a  stream  function,  the  conservation  of  mass  is  satis¬ 
fied  identically.  That  is,  'I'  is  defined  such  that 
d* 


3x 


=  -prv 


(5.5) 


—  ■  »v 


If  we  take in  the  form 


*  <£,  n)  =^2 7s  f 


(5.6) 


(5.7) 


and  introduce  the  Mangler-Dorodnitsyn  transformation,  Refs.  (44),  (45)  and  (46): 


T)  = 


P  u 
e  e 


^2~T 


j  r  ^ ■ -  dy 
J  o  p 


\  -  o  p 

o  e 


JA  2 

p  n  u  r  dx 
w  w  e  o 
o 


(5.8) 


(5.9) 


the  conservation  equations  are  transformed  into  a  set  of  non-linear  ordinary  differential 
equations  in  the  following  form:  The  conservation  of  momentum  becomes 
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The  conservation  of  species  takes  on  the  form 
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The  conservation  of  energy  for  an  optically  thick  gas  becomes 
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It  is  noted  that  the  following  definitions  were  utilized  in  deriving  the  preceding  equations. 
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The  boundary  conditions  to  be  utilized  with  Eqs.  (5. 10),  (5. 11)  and  (5. 12)  are  as  follows. 
At  the  surface,  assuming  zero  slip  and  no  mass  transfer: 
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(5. 14) 


(5. 15) 


(5. 16) 


The  surface  boundary  conditions  on  chemical  composition  are  given,  in  general,  by 
appropriate  chemical  kinetic  equations.  If  the  gas  is  assumed  to  be  in  equilibrium  at 
the  surface,  one  simply  utilizes  the  appropriate  equilibrium  constants. 

At  the  outer  edge  of  the  boundary  layer,  one  may  apply  the  asymptotic 


30 


I 


I 


conditions: 


lim  f  =  lim  0  =  1. 0 
T)  - 00  ^  rj  -® 


(5. 17) 


and  in  addition,  one  sets  the  following  boundary  condition  on  the  chemical  composition: 
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5. 2  Optically  Thin  Boundary  Layer 


(5. 18) 


When  the  boundary  layer  is  optically  thin,  one  may  utilize  the  equations 

for  the  conservation  of  mass,  momentum,  and  chemical  species  given  by  Eqs.  (5.10) 

and  (5. 11),  however,  utilizing  Eqs.  (3.  32)  and  (3. 35)  the  conservation  of  energy  now 
takes  on  the  form: 
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which  may  be  transformed  into  the  form: 
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Although  the  boundary  conditions  at  the  surface  are  the  same  as  those  given  in  Eqs. 

(5. 14),  (5. 15)  and  (5. 16),  the  edge  boundary  conditions  are  applied  somewhat  differently 
in  that  asymptotic  behavior  is  not  necessarily  expected.  Instead,  one  has  a  floating 
boundary  condition,  in  much  the  same  fashion  as  the  low  Reynolds  number  layer 
problem,  Ref.  (47)  where  an  auxiliary  integral  constraint  is  introduced.  That  is, 
here  one  may  write  in  first  approximation: 

6(5)  =  f  (6)  =  1.0;  C.(6)  =  C.  (5.21) 

e 

where  6,  the  thickness  of  the  radiating  viscous  layer  is  not  known  "a  priori"  but  may 
be  determined  by  means  of  additional  constraints  which  are  based  on  the  requirement 
that  each  variable  and  its  derivatives  should  match  smoothly  with  the  outer  shock 
layer  solution. 

5. 3  Gas  Properties 

Two  different  gas  models  were  treated  in  the  numerical  computations,  constant 
properties  (Figures  12,  13)  and  variable  properties  (Figures  14  and  15).  For  the  former, 
it  was  assumed  that  the  product  of  viscosity  and  density,  as  well  as  the  mean  molecular 
weight  was  a  constant  across  the  boundary  layer;  for  the  latter,  these  were  calculated 
locally  and  hence  were  variable,  (see  Ref.  49  and  TABLE).  The  Prandtl  and  Lewis  numbers 
were  assumed  to  be  equal  to  unity  in  both  cases,  and  for  simplicity  the  thermal  conducti¬ 
vity  was  assumed  to  be  proportional  to  the  square  root  of  the  absolute  temperature. 


32 


6.  DISCUSSION  OF  RESULTS 


In  this  section,  we  will  discuss  the  numerical  results  obtained  during 
the  course  of  the  present  investigation.  We  begin  by  remarking  that  the  motivation 
in  carrying  out  the  study  was  to  develop  additional  understanding  of  radiative  transport 
processes  in  the  hypersonic  shock  layer.  Referring  to  Figure  1,  we  observe  that  the 
shock  layer  has  been  characterized  here  as  consisting  of  three  distinct  regions: 

Region  I  -  Shock  wave 

5k 

Region  II  -  Influscid  layer 

Region  HI  -  Boundary  layer 

Since  the  equilibration  of  energy  is  accomplished  throughout  the  shock 
layer,  the  "temperature"  has  been  represented  initially  by  four  separate  curves  re¬ 
presenting  the  translational,  rotational,  vibrational  and  electronic  excitational  modes. 
After  thermal  equilibration  is  complete,  of  course,  the  four  separate  curves  merge 
into  a  single  curve  of  temperature.  It  is  seen  that  the  normal  component  of  velocity 
decreases  monotonically  as  the  flow  approaches  the  vehicle,  and  in  the  absence  of  abla¬ 
tion  or  other  forms  of  mass  transfer,  vanishes  identically  at  the  surface.  It  is  also 
seen  that  the  gas  density  rises  as  the  flow  approaches  the  surface.  Simultaneously 
with  the  change  in  the  thermodynamic  properties,  a  change  in  chemical  composition 
is  also  taking  place.  In  fact,  even  before  thermal  equilibration  is  necessarily  com¬ 
pleted,  the  chemical  processes  of  dissociation  and  ionization  begin.  But  this  clearly 

This  word  represents  the  first  author's  generalization  of  the  word  "inviscid"  to 
include  the  absence  of  viscous  effects,  diffusion  effects  and  thermal  conduction. 
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depends  on  the  intermolecular  collision  frequency  and  the  rates  of  the  various  chemi¬ 
cal  kinetic  processes.  In  general,  the  problem  is  extremely  complicated  since  (me 
is  dealing  with  a  number  of  complex  physicochemical  processes  some  of  which  take 
place  in  parallel  and  others  in  series.  Finally,  when  one  introduces  the  additional 
complication  of  radiation,  which  is  in  turn  related  to  the  gas  temperature,  the  number 
density  of  chemically  reacting  species,  and  the  population  of  states  involved  in  radiative 
transitions,  the  total  problem  almost  appears  to  be  hopeless.  Nevertheless,  great 
progress  is  being  made,  and  in  this  paper  it  is  expected  that  additional  understanding 
of  coupled  radiative  phenomena  will  be  added  to  the  existing  store  of  information. 

It  is  now  well  known  that  the  physical  extent  of  each  of  these  layers 
is  a  function  of  the  Reynolds  number,  Mach  number,  Knudsen  number,  surface  tem¬ 
perature  and  vehicle  geometry.  Thus,  at  low  altitudes.  Regions  I  and  III  are  infinite¬ 
simally  thick,  and  Region  II  envelopes  the  vehicle,  whereas  during  high  altitude  flight, 
the  two  viscous  Regions  (I  and  III)  expand,  merge  and  "swallow"  Region  II.  In  the 
present  study,  we  have  focussed  our  attention  on  Regions  I  and  in,  for  several  reasons, 
including  the  fact  that  these  regions  are  least  well  understood,  and  further,  where 


Region  n  dominates  the  flow,  the  existing  theoretical  models  seem  quite  reasonable. 

Figures  2  and  3  show  the  values  used  by  us  for  the  Rosseland  mean  absorp¬ 
tion  coefficient  for  treating  an  air-like  gas  in  the  optically  thick  limit.  These 
may  be  represented  analytically  by  the  following  expressions  for  the  linearized  and 
quadratic  representations: 
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Kj,  =  4.52  X  10-7  (p)1'31  4<6-18  *  10'4t  -  7-13  *  «'V)  (6>„ 

where  the  Rosseland  mean  absorption  coefficient  1^  is  expressed  in  cm. -1,  the  pressure 

p  in  atmospheres,  and  the  temperature  T  in  °K.  In  our  computations  for  an  optically 

thin  gas,  we  used  the  same  functional  form  for  the  Planck  mean  absorption  coefficient, 

but  simply  assumed  that  K  =  8.  3K  . 

P  R 

In  Figures  4,  5  and  6  are  depicted  solutions  obtained  for  the  structure 
of  a  normal  shock  wave  in  an  optically  thick  radiating  gas.  The  range  of  Mach  numbers 
treated  was  12. 5  <  M  <  50.  The  symbol  a^  denotes  the  maximum  value  of  the  radia¬ 
tion  parameter  discussed  in  Appendix  A.,  that  is;  using  symbols  of  Eq.  (A. 9): 
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In  this  equation,  KR  is  given  by  Eq.  (6. 1),  T  is  given  by  Eq.  (A.  10)  and  aj  is 
shown  in  Figure  19  as  a  function  of  P,  defined  in  Eq.  (4. 4).  With  a.^  =  0,  no 
radiation  is  present.  As  a^  is  arbitrarily  increased,  it  is  seen  that  radiation  broad¬ 
ening  occurs  on  the  upstream  side  of  the  shock  wave.  These  results  are  qualita¬ 
tively  similar  to  those  obtained  by  Sen  and  Guess,  Ref.  (28)  who  treated  the  low 
speed  range  of  Mach  numbers  1.  5  <  M  <  4,  in  that  the  radiation  broadening 

I 

occurred  in  the  low  temperature  (upstream)  side  of  the  shock  wave.  However, 
since  their  choice  of  transport  and  radiative  properties,  as  well  as  Mach  number 
range,  was  different  (see  Appendix  A),  their  quantitative  results  are,  of  course, 
different. 

Upon  introducing  the  following  definition  of  the  slope  thickness  in  the 
absence  of  radiation, 
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we  obtain  a  reference  value  of  the  nominal  shock  thickness.  In  order  to  calculate 
the  nominal  increase  in  shock  wave  thickness  as  a  function  of  Mach  number  and  radia¬ 
tion  function  a^,  the  following  arbitrary  procedure  was  adopted.  First,  we  write: 


(AY)ai  =  (Y2-Yi)ai 


The  point  (  Y  )  is  defined  as  the  point  of  intersection  of  the  slope 

1  u 


[(-£-) ...  1 


and  the  asymptote  T  .  We  make  a  vertical  projection 


ai=° 


from  the  point  (Y1>a  _  0  up  to  the  curve  T  =  KL  _  0  which  defines  a  refereni 
temperature  at  the  point  (Y  )  .  For  each  solution  curve,  the  point  (Y  )  , 

J.  U  X  Ej f  U 

was  determined  at  this  reference  temperature.  Since  the  results  on  the  downstream 

side  of  the  shock  wave  are  insensitive  to  a  ,  Y  is  defined  as  the  fixed  point  given  by: 
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which  defines  a  reference 
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which  is  also  the  intersection  of  the  slope 
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with  the  asymptote 
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The  ratio  AY/  A Yq  is  shown  in  Figure  7,  where  it  is  seen  that  radiation 
broadening  is  less  pronounced  at  higher  Mach  numbers  for  a  given  value  of  a^  . 


It  is  noted  that  in  this  model  of  shock  wave  structure  in  an  optically  thick 


I 


gas,  we  have  assumed  that  the  gas  is  characterized  by  a  single  temperature;  however, 
it  has  not  been  assumed  that  local  thermochemical  equilibrium  exists,  in  that  the 
chemical  composition  was  assumed  to  be  frozen  at  the  upstream  value.  No  attempt 
was  made  to  include  chemical  kinetics,  and  thus,  the  results  would  tend  to  apply  to 
the  initial  part  of  a  shock  layer  which  is  optically  thick.  In  order  to  be  self-consistent, 
we  should  have  used  the  absorption  coefficient  Corresponding  to  this  frozen  chemical 
composition,  (e.g. ,  in  air,  corresponding  to  the  molecular  bands  of  N  and  O  ), 

Z  La 

rather  than  that  corresponding  to  complete  thermochemical  equilibrium  (approximately 
given  in  Figures  2  and  3).  However,  estimates  indicated  that  the  difference  between 
the  two  was  not  too  large  for  most  of  the  conditions  of  interest  here.  Thus,  in  this  para¬ 
metric  study,  we  used  the  values  of  given  in  Figure  2. 

The  numerical  results  obtained  for  the  other  extreme  limiting  case  of 
shock  wave  structure,  in  an  optically  thin  radiating  gas,  are  shown  in  Figures  8  through 
11.  Two  values  of  the  Mach  number,  M  =  25  and  35  were  treated.  Now,  as  noted, 
in  an  optically  thin  layer,  the  photon  mean  free  path  is  large  compared  with  the  thickness 
of  the  layer,  and  hence,  one  may  assume  that  the  gas  emits  but  does  not  absorb  radia¬ 
tive  energy.  Consequently,  unlike  the  optically  thick  shock  wave  where  the  radiative 
effects  appear  on  the  upstream  (low  temperature)  side,  here,  the  radiative  effects 
appear  primarily  on  the  downstream  (high  temperature)  side  of  the  shock  wave.  That 
is,  the  high  temperature  region  of  the  gas  radiates  in  both  the  upstream  and  downstream 
directions,  but  since  energy  escapes  from  the  structure  of  the  shock  wave,  the  peak 
temperature  which  is  reached  is  lower  than  would  have  been  obtained  if  radiation  effects 
were  absent. 
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In  this  part  of  the  study,  the  symbol  a^  was  introduced  as  a  measure  of 
the  magnitude  of  the  radiative  term  appearing  in  the  equation  for  the  conservation  of 
energy  in  an  optically  thin  shock  wave  (see  Section  4).  Thus,  the  case  a  =  0  corres- 
ponds  to  zero  radiation,  while  the  effects  of  radiation  are  enhanced  as  a0  is  increased. 
Examination  of  Figure  8  shows  that  as  a  is  increased,  two  things  happen;  one  is  that 
the  peak  temperature  is  decreased,  but  moreover,  the  downstream  conditions  are  not 
fixed  constants  determined  "a  priori"  from  the  Rankine-Hugoniot  conditions.  Further, 
the  downstream  conditions  are  not  necessarily  approached  asymptotically,  but  the  tem¬ 
perature  continues  to  decay  with  increasing  distance.  That  is,  since  the  high  tempera¬ 
ture  gas  is  continually  radiating  energy  away  as  it  flows,  the  temperature  continues  to 
drop,  and  the  solution  must  be  matched  to  that  obtained  in  Region  II.  Of  course,  if 
Region  n  does  not  exist,  due  to  the  viscous  merging  of  Regions  I  and  III,  a  matching 
procedure  must  be  utilized  at  the  interface  between  the  contiguous  layers. 

When  the  effects  of  dissociation  equilibrium  are  included,  the  chemical 


energy  required  for  dissociation  acts  to  decrease  the  temperature  throughout  the  shock 
wave.  Hence,  the  radiative  transport,  whose  temperature  dependence  appears  explictly 

4 

as  T  and  implicitly  through  K  =  K  (T,  p),  is  much  less  important.  Consequently, 

P  P 

at  a  Mach  number  of  25,  and  for  the  same  values  of  a  as  utilized  in  the  non-dissociating 
model,  i. e. ,  0^  a  £  10,  it  is  seen  (see  Figure  9)  that  there  are  no  significant  coupled 
radiative  effects  for  a  ^  10. 

On  the  other  hand,  at  the  higher  Mach  number  (M  =  35),  even  when  disso¬ 
ciation  equilibrium  was  included  in  the  model,  (see  Figures  10  and  11),  radiative 


coupling  was  significant  for  a  value  of  a  =  10.  The  specific  equation  for  a  is: 

4  2 
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where  the  symbols  are  introduced  in  Appendices  A  and  B.  It  was  found  that  the 
solution  was  relatively  insensitive  to  the  effect  of  pressure  on  dissociation  and 
hence  the  error  resulting  from  the  fact  that  a  particular  choice  of  a^  also  specifies 

i 

a  pressure  is  small. 

Here  again,  it  is  noted  that  a  simplified  model  of  the  chemistry  was 
utilized  in  order  to  obtain  qualitative  results  in  a  simple  manner.  Thus,  since  it 
was  arbitrarily  assumed  that  the  dissociation  process  was  rapid,  and  the  ionization 
process  was  very  slow,  the  results  are  applicable  to  the  initial  portion  of  an  optically 
thin  shock  wave.  Since  the  ionization  process  also  acts  to  reduce  the  gas  temperature 
by  absorbing  the  available  energy,  the  radiative  coupling  might  not  actually  be  as  pro¬ 
nounced  as  shown  in  Figures  10  and  11.  It  should  perhaps  also  be  mentioned  that  as  in  the 
case  of  the  optically  thick  shock,  and  as  done  in  the  treatment  of  the  boundary  layer 
as  well,  we  used  the  functional  form  of  the  absorption  coefficient  corresponding  to  local 
thermochemical  equilibrium.  That  is,  we  did  not  take  into  account  the  effect  of  de¬ 
partures  from  chemical  equilibrium  on  the  functional  form  of  and  K^. 

In  Figures  12  through  16  are  presented  the  numerical  results  of  an  investi¬ 
gation  of  the  optically  thick  boundary  layer,  (Region  HI),  for  two  values  of  the  boundary 
layer  edge  temperature,  Tg  =  20,000°R  and  40,000°R.  The  wall  temperature  was  arbi¬ 
trarily  taken  equal  to  4, 000°R  and  the  computations  were  carried  out  for  the  case  of 
constant  properties  and  variable  properties.  Curves  of  the  variation  of  tangential  velo¬ 
city  u  and  gas  temperature  T,  normalized  by  their  edge  values,  are  shown  as  a  function 
of  tj,  the  stretched  y  coordinate.  In  this  part  of  the  study,  the  symbol  ag  denotes  the 
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maximum  value  of  the  radiation  parameter,  which  therefore  plays  the  same  role  as 
a  ,  where  ag  is  given  by  Eq.  (A.  4)  and  T  is  obtained  from  Eqs.  (A.  6)  and  (A.  7).  The 
dependence  of  a  upon  p  is  shown  in  Figure  19  for  both  the  linear  and  quadratic 

u  G 

forms  of  K  . 


Here  again,  it  is  seen  that  radiation  broadening  occurs  as  a  is  increased. 

O 

However,  as  shown  in  Figure  16,  thickening  of  the  layer  is  not  the  only  effect.  Upon 
making  a  plot  of  the  ratio  of  the  sum  of  the  conductive  and  radiative  terms  normalized 
by  the  conductive  term,  where  the  numerical  solutions  were  obtained  by  keeping 
the  boundary  conditions  fixed  and  merely  changing  a  ,  one  finds  that  this  ratio  in- 

O 

creases  linearly  with  a  ,  and  at  a  rather  substantial  rate. 

In  Figures  17  and  18,  we  have  shown  a  typical  set  of  numerical  results 
obtained  for  an  optically  thin  boundary  layer.  The  variation  of  the  dimensionless 
tangential  velocity  and  gas  temperature  is  shown  as  a  function  of  the  variable  77, 
for  three  different  values  of  the  radiation  parameter  a  .  Here  it  is  noted  that  a  is 


given  by: 
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which  may  be  simplified  by  introducing  the  hypersonic  approximation 
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so  that  we  also  obtain: 
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(6.9) 


(6. 10) 


(6.11) 


Which  one  of  the  solutions  given  in  Figures  17  and  18  applies  in  a  given  physical 
situation  would  depend  on  matching  this  solution  smoothly  into  the  solution  to  the 
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outer  contiguous  layer. 

The  condition,  Eq.  (3. 14) ,  for  the  shock  layer  to  be  optically  thin  can 

always  be  reached  for  sufficiently  small  values  of  the  vehicle  nose  radius  R  .  The 

B 

values  pertaining  to  our  choices  for  a  and  a  ,  given  by  Eqs.  (6. 7)  and  (6. 11)  are, 

£*  4 

for  the  most  part,  physically  reasonable.  On  the  other  hand,  the  pressures  corres¬ 
ponding  to  the  choices  for  a  and  a  are  seen  from  Figure  19  to  be  extremely  large, 

1  i) 

5  6 

p  ~  10  to  10  atmospheres.  Moreover,  condition  (3. 37)  (which  as  noted  in 
Section  3  implies  that  condition  (3. 13)  applies  as  well  for  the  shock  wave  or  boundary 
layer)  is  found  by  inspection  of  our  results  not  to  be  satisfied  when  KR  is  intro¬ 
duced  for  ,  This  indicates  that  the  shock  wave  and  boundary  layer  are  not  optically 
thick.  In  fact,  the  total  shock  layer  is  usually  more  nearly  optically  thin  than  thick. 
However,  upon  introducing  the  hypersonic  approximation  that  the  total  shock  layer 
thickness  is  of  the  order  of  0. 1  R  ,  we  see  from  Figure  3  that  condition  (3. 13)  can 
be  approximately  satisfied  for  the  total  shock  layer  for  realistic  values  of  the  stagnation 
pressure  and  large  values  of  R  ,  such  as  may  be  expected  for  manned  space 
vehicles.  If  the  total  shock  layer  is  optically  thick  and  the  gas  temperatures  large 
enough  for  radiation  to  dominate  over  aerodynamic  transport  processes,  one  would 
expect  the  boundary  layer  and  shock  wave  to  adjust  themselves  such  that  they  would 
each  be  of  the  order  of  a  photon  mean  free  path  in  thickness,  i.e. , 


j^Kdy-l  (6.12) 

where  6  is  the  thickness  of  the  shock  wave  or  boundary  layer  and  K  an  average  of  K^over  v. 
The  reason  for  expecting  this  self-adjustment  in  the  thickness  is  that  if  the  shock 
wave  and  boundary  layer  were  thin  while  the  total  layer  is  thick,  black  body  radiation 
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flux  would  flow  out  each  side  of  the  shock  layer.  If  the  temperature  is  high, 
the  heat  content  of  the  shock  layer  would  not  be  great  enough  to  sustain  this  and 
the  shock  wave  and  boundary  layer  would  broaden  out  to  the  extent  that  a  signifi¬ 
cant  portion  of  the  radiative  flux  from  Region  n  would  be  reabsorbed. 


:  i 

i  1 

t  i 

i 


7.  CONCLUSIONS 

In  this  paper,  which  gives  an  introductory  treatment  of  the  coupling  between 
radiative  transport  processes  and  hypersonic  chemically  reacting  flows,  four 
different  problems  have  been  investigated,  i.  e.  the  hypersonic  shock  wave  structure 
problem  and  the  hypersonic  laminar  boundary  layer  problem,  for  both  the  optically 
thin  and  optically  thick  limits. 

The  true  physical  situation  will  lie  between  these  limits,  and  usually  much 
closer  to  the  optically  thin  limit,  however,  by  treating  the  problem  in  the  two  limits,  it  is 
felt  that  certain  features  of  the  true  solution  have  been  captured.  For  example, 
in  the  case  of  shock  wave  structure,  it  is  anticipated  that  the  peak  shock  temperature 
will  be  diminished  by  radiation,  since  some  of  the  radiation  escapes  without  being 
reabsorbed.  On  the  other  hand,  some  of  the  radiation  is  absorbed,  and  hence,  the 
shock  wave  will  be  broadened  on  the  upstream  side. 
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TABLE  -  GAS  PROPERTIES 


c  -o.a-BHL 

p  lb  °R 


hN  =  1. 43  x  104  +  0.  35  T(°R), 


h =-280  +  0.  32  T(°R), 
2 


k  =  2.  79  x  10“5  Y  T(°R)  ,  ■■  Watt8 

cm  °K 


a  =  5. 67  x  10 


-12  watts 

2  o4 
cm  K 


. -5  /  717  \  /  T(°H)  \  3/2  lb 

“',-16x10  W  tsr)  ‘srsr 

Note  that  k  was  utilized  in  the  calculation  of  the  radiation  parameter  and  that  /j  was 
utilized  in  calculating  l ,  but  not  in  the  calculation  of  the  Prandtl  number  which  was 
assumed  to  be  equal  to  unity. 
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APPENDIX  A  -  COMPARISON  OF  RADIATION  AND  CONDUCTION  TERMS  IN  AN 
OPTICALLY  THICK  GAS 

For  an  optically  thick  gas,  we  may  write  the  energy  flux  vector  in  the 

form: 


S=-k 


16  C7T3\ 
3kKR  ) 


VT  +  L  Pl  Vtht 


(A.  1) 


In  order  to  assess  the  relative  importance  of  the  radiative  and  con¬ 
ductive  contributions  in  both  a  boundary  layer  and  a  shock  wave,  we  must  compare 

3 

the  magnitude  of  the  dimensionless  radiative  function  16  O  T  /3kK  with  unity.  Since 

K 

the  former  is  not  a  constant,  the  comparison  should  be  made  at  that  value  of  the  tem¬ 
perature  at  which  the  radiative  function  is  a  maximum. 

Let  us  therefore  introduce  the  following  functional  forms: 


k  =  aT 


and 


b  cT 


Kn  =  |3p  e 


-dT2 


(A.  2) 


(A.  3) 


! 

I 

I 

f 

i 


1 

i 


I 

\ 

\ 
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where  a  ,  a,  b,  c  and  d  are  constants.  Then,  there  follows: 


16  aT3 
3kKR 


16  a 


3aj8  p 


TmedT 


CT 


(A.  4) 


where  m  =  3  -  a.  If  c  >  >  d,  the  radiative  function  will  be  a  maximum  at  the  lower 
gas  temperatures  which  appear  in  the  problem,  and  will  decrease  as  the  temperature 
rises.  In  general,  however,  the  radiation  function  will  not  have  a  monotonic  behavior, 
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and  therefore  we  are  interested  in  determining  the  temperature  at  which  a  maximum 
is  attained. 

In  a  boundary  layer,  the  local  pressure  is  a  constant,  hence  upon 
setting  the  partial  derivative  of  Eq.  (A.  4)  with  respect  to  temperature,  at  constant 
pressure,  equal  to  zero,  we  obtain  the  quadratic  equation 

2dT2-cT  +  m  =  0  (A.  5) 


and  hence  the  radiative  function  is  stationary  when  the  temperature  is  given  by 
T  =  m/c  ;  (d  =  0) 


T  = 


8dm 


(d/  0) 


(A.  6) 
(A.  7) 


2 

If  m  >  0,  and  c  >  8  dm,  where  a,  b,  c  and  d  are  positive  constants,  it  is  seen  that 
the  radiation  function  goes  through  a  maximum  when  one  utilizes  the  minus  sign  (low 
gas  temperatures)  and  goes  through  a  minimum  when  one  selects  the  plus  sign  (high 
gas  temperatures). 

In  the  case  of  the  shock  wave  structure  problem,  we  note  that  the 
flow  is  one-dimensional,  and  moreover,  the  pressure  is  not  constant  through  the 
layer.  Hence  we  may  introduce  the  equation  of  state  and  the  conservation  of  mass 
to  obtain 


in  which  case  the  radiation  function  becomes: 


160  T3 
3kKR 


(A.  9) 
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where  n  =  3  -  a  -  b.  If  one  assumes  that  both  It  and  v  are  changing  slowly  in  the  low 
temperature  regime,  then  the  radiative  function  is  a  maximum  when 

T  «  n/c;  d  =  0  (A.  10) 

T  ■"  8-<jD  X  ;  (d  i-  0)  (A.  11) 

Thus,  upon  substituting  either  Eq.  (A.  6)  or  (A.  7)  into  Eq.  (A.  4),  or 
either  Eq.  (A.  10)  or  (A.  11)  into  Eq.  (A.  9),  respectively  one  can  immediately  determine, 

"a  priori",  the  relative  importance  of  the  radiation  and  conduction  terms.  It  should 

be  noted  that  if  the  radiative  function  is  a  sharply  peaked  function  of  temperature,  then, 

of  course,  the  average  value  of  the  radiative  function  can  be  considerably  smaller  than 
its  maximum  value. 

It  is  also  noted  that  if  we  had  selected ,  for  example,  the  functional 
forms  for  thermal  conductivity  and  absorption  coefficient  utilized  by  Sen  and  Guess 


(Ref.  28)  i.  e. 

k  =  a  T2'5 

(A.  12) 

kr  =  p  pt"3'5 

(A.  13) 

3 

then  the  radiative  function  would  increase  monotonically  very  nearly  as  T  and  the 
major  contribution  of  the  radiation  would  appear  at  the  highest  temperatures  en¬ 


countered  in  a  given  physical  problem. 
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APPENDIX  B  -  NORMALIZATION  OF  THE  SHOCK  WAVE  EQUATIONS 


We  will  derive  here  the  non-dimensional  form  of  the  equations.  Following 
Grad's  notation,  Ref.  (48),  we  introduce  the  following  dimensionless  variables 


mv  Po 

1  \  ~  — — —  •  V  =  ~~  ■■■ 

P  ’  2 

m 

2_  ^ 
m  R  T 

m  _  p 

T  p2  5  CT  P 


5=  “S- 

P2 


Upon  defining  the  transport  lengths. 


Lt  = 

m  t 


and  putting  the  transport  properties  in  the  form 


M  =  CjT 

#  =  c  t3/^ 

AM  2  /v 


k=  c3t 


We  can  select  a  reference  temperature  T  and  introduce  the  non-dimensional  spatial 
variable 


where 
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v 


2  * 
m  R  T 
*  m 

T  =  - — — 


(B-5) 


and 


* 

L 


m 


(B-6) 


The  following  equations  are  then  obtained  in  terms  of  the  non-dimensional  variables. 
The  conservation  of  mass  is: 


ru  =  1.0 


The  conservation  of  momentum  is: 


Tt1/2u  u2+  T(1  +  Ca)-u 


The  conservation  of  species  is: 


*AU*A-*<CA-CA  > 


where  the  diffusion  flux  is: 


,  ^  r_j 

d  L  i  + 


dC, 


V*  -  °A> 


C.  dY 
A 


1  +  C  ■.  K 
A 


T  dr  1 

T  dY  J 


and  the  chemical  source  term  is: 

* 

L 


= 


A  *1/2 
T 


dY 


1 

5 


do 

dY 


For  an  optically  thick  gas,  the  conservation  of  energy  is  given  by: 
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(B-7) 


(B-8) 


( B-9) 


(B-10) 


(B-ll) 


f  I 


1  + 


16  CT 


3kK 


1/2 


R  J 


dT  J:A  [/-  -  \  A„  kT(1 

"dlT  =  "m"  [  °A  ”  Cm'T  +  *H  + 


+  CA>  ] 

i-cA)  TJ 


CA*  -A 


+  C  , 


L 
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In  the  above  equations: 
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and  the  dimensionless  ratios  are 
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where  0  has  the  behavior  of  the  reciprocal  of  the  Prandtl  number  and  0.  has  the 

t  pCr 

behavior  of  the  reciprocal  of  the  Schmidt  number. 

In  addition  one  has  the  equation  of  state  for  a  dissociating  diatomic  gas 
cr  =  r  (1  +  C  )  T  (B-15) 

For  an  optically  thin  radiating  gas,  equation  (B-12)  is  replaced  by: 
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where  the  dimensionless  absorption  coefficient  is  given  by: 
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(B-17) 


57 


SCHEMATIC  REPRESENTATION  OF 
HYPERSONIC  SHOCK  LAYER 


FIGURE  I 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THICK 

RADIATING  GAS,  M  =  I2.5 


FIGURE  4 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THICK 

RADIATING  GAS,  M  =  25 
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FIGURE  5 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THICK 

RADIATING  GAS,  M  =  50 
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FIGURE  6 


DEPENDENCE  OF  NOMINAL  SHOCK  WAVE  THICKNESS 


FIGURE 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THIN 

RADIATING  GAS,  M  =  25 


FIGURE 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THIN 


FIGURE 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THIN 

RADIATING  GAS,  M  =  35 


FIGURE  10 


SHOCK  WAVE  STRUCTURE  IN  AN  OPTICALLY  THIN 
RADIATING  GAS  WITH  DISSOCIATION,  M  =  35 
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FIGURE  II 


BOUNDARY  LAYER  STRUCTURE  IN  AN  OPTICALLY  THICK 
RADIATING  GAS,  (CONSTANT  PROPERTIES),  Te  =  20, 000°R 


FIGURE  12 


BOUNDARY  LAYER  STRUCTURE  IN  AN  OPTICALLY  THICK 
RADIATING  GAS,  (CONSTANT  PROPERTIES),^  =40,000°  R 


FIGURE  13 


BOUNDARY  LAYER  STRUCTURE  IN  AN  OPTICALLY  THICK 
RADIATING  GAS,  (VAR  I  ABLE  PROPERT I  ES),Te  =  20,000°  R 


BOUNDARY  LAYER  STRUCTURE  IN  AN  OPTICALLY  THICK 


r 


FIGURE  15 


NORMALIZED  RADIATIVE  PLUS  CONDUCTIVE 

HEAT  TRANSFER 


FIGURE  16 


VARIATION  OF  RADIATION  PARAMETER  FOR 


GENERAL  Hi  ELECTRIC 
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The  effect  of  radiative  transport  on  the  structure  of  the 
Bhock  layer  produced  In  roni  of  a  hypersonic  vehicle  in  an  earth- 
,ike  atmosphere  is  studied.  Specifically,  the  effects  of  radiation 
on  the  structure  of  the  boundary  layer  and  the  shock  wave  are 
determined  numerically  :or  the  optically  thick  and  optically  thin 
Limits, 


In  tiW  optically  thin  limit,  the  effect  of  radiation  is  to  reduce 
the  peak  shock  tempertiure.  Also,  on  the  downstream  side,  the 
temperature  continues  tc  decay  with  distance  because  the  gas 
continues  to  emit,  and  tte  solution  must  be  matched  to  that  per¬ 
taining  to  the  contiguous  inner  viscous  layer.  On  the  other  hand, 
the  peak  shock  temperature  and  the  Rankine-Hugoniot  conditions 
are  unaffected  by  the  radiation  in  the  optically  thick  limit,  but 
both  boundary  layer  and  shock  wave  are  broadened  by  the  radia¬ 
tion.  Although  it  is  found  that  the  boundary  layer  and  shock  wave 
are  n«»h««Hy  never  optically  thick,  neither  are  they  always  opti¬ 
cally  Thus,  some  of  the  radiation  emitted  is  reabsorbed  in 

die  boundary  layer  and  shock  wave  and  they  are  thereby  "radia¬ 
tion"  broadened  to  some  degree.  The  broadening  is  expected  to 
be  largest  for  high  temperatures  and  when  the  total  shock  layer 
is  opticslly  thick  or  nearly  so.  It  is  also  pointed  ™it  that  the  gi ay 
gas  approximation  is  very  poor  for  high  temperature  air. 
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